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Prognostics technologies determine the health (or damage) state of a component or sub- 
system, and make end of life (EOL) and remaining useful life (RUL) predictions. Such infor- 
mation enables system operators to make informed maintenance decisions and streamline 
operational and mission-level activities. We develop a model-based prognostics method- 
ology for pneumatic valves used in ground support equipment for cryogenic propellant 
loading operations. These valves are used to control the flow of propellant, so failures may 
have a significant impact on launch availability. Therefore, correctly predicting when valves 
will fail enables timely maintenance that avoids launch delays and aborts. The approach 
utilizes mathematical models describing the underlying physics of valve degradation, and, 
employing the particle filtering algorithm for joint state-parameter estimation, determines 
the health state of the valve and the rate of damage progression, from which EOL and 
RUL predictions are made. We develop a prototype user interface for valve prognostics, 
and demonstrate the prognostics approach using historical pneumatic valve data from the 
Space Shuttle refueling system. 


Nomenclature 


ISHM 

Integrated Systems Health Management 

EOL 

end of life 

RUL 

remaining useful life 

LH2 

liquid hydrogen 

t or k 

time 

tp or kp 

time of prediction 

X 

state vector 

e 

parameter vector 

V 

process noise vector 

y 

output vector 

n 

measurement noise vector 

Teol 

EOL threshold 

RMAD 

relative median absolute deviation 

RA 

relative accuracy 


I. Introduction 

Within Integrated Systems Health Management (ISHM), technical advances in and increased acceptance 
of automated diagnosis technologies are leading to a shift in research focus to the area of automated prognos- 
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tics. Prognostics algorithms focus on determining the health (or damage) state of a component or subsystem, 
and making end of life (EOL) and remaining usefid life (RUL) predictions. EOL and R.UL predictions en- 
able system operators to make informed maintenance decisions and streamline operational and mission-level 
activities. 

As a case study, we consider ground support equipment from the Space Shuttle cryogenic refueling system. 
This system transfers cryogenic propellants from a storage tank to the vehicle tank through a network of 
pipes and valves. Since the valves are used to control the flow of propellant, failures may have a significant 
impact on launch availability. Hence, valve prognostics can provide significant value to launch operations. 
We focus on one of the pneumatic transfer line valves in this system, and apply a model-based prognostics 
approach. 

For the selected valve, we develop a complete prognostics solution, including models, algorithms, and a 
user interface. The prognostics algorithms are model-based, meaning that, unlike data-driven approaches, 1 
we use mathematical models describing the underlying physical principles of valve operation and how valves 
fail. 2-6 We construct a detailed physics-based model of a pneumatic valve that includes models of different 
damage mechanisms and how they progress in time. Using this model, we employ the particle filtering 
algorithm for joint state-parameter estimation to determine the health state of the valve and the rates of 
damage progression, from which EOL and RUL predictions are made. Although we focus on a pneumatic 
valve, the approach is general in that it may be applied to any component, given its model. We also develop 
a prototype user interface to communicate the algorithm results, based on an analysis of what information 
the algorithms provide and how the information should be made available to the user. We demonstrate the 
prognostics approach as applied to a real-world case study using historical pneumatic valve data from the 
Space Shuttle refueling system. 

The paper is organized as follows. Section II formulates the prognostics approach. Section III describes 
the valve model. Section IV describes the damage estimation algorithm, and Section V describes the pre- 
diction algorithm. Section VI develops the prognostics interface. Section VII demonstrates the approach on 
historical pneumatic valve data. Section VIII concludes the paper. 


II. Prognostics Approach 

The goal of prognostics is to predict the EOL and/or the RUL of a component. In this section, we first 
formally define the problem of model-based prognostics. We then describe a general model-based architecture 
within which a prognostics solution may be implemented. 


II. A. Problem Formulation 

We assume the system may be described by 

x(f) = f(f,x(t),0(f),u(t),v(f)) 

y(t) = h(t,x(f),0(f),u(t),n(t)), 


where x(f) £ R”* is the state vector, 9(f) £ R" 0 is the parameter vector, u(<) £ R" u is the input vector, 
v(t) £ R"" is the process noise vector, f is the state equation, y (t) £ R" B is the output vector, n(f) £ R™’ 1 
is the measurement noise vector, and h is the output equation. This form represents a general nonlinear 
model with no restrictions on the functional forms of f or h. Further, the noise terms may be coupled in 
a nonlinear way with the states and parameters. The parameters 9(f) evolve in an unknown way, but, in 
practice, are typically considered to be constant. 

Our goal is to predict EOL at a given time point tp using the discrete sequence of observations up to 
time tp , denoted as yo-.t P - EOL is defined as the time point at which the component no longer meets a 
functional requirement (e.g., a valve does not open in the required amount of time). This point is expressed 
through a damage threshold, beyond which the component fails to function properly. In general, we may 
express this threshold as a function of the system state and parameters, T E OL(x(t), 9(f)), which determines 
whether EOL has been reached, where 


T E OL(.x(t),9(t)) 


1, if EOL is reached 
0, otherwise. 
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Using this function, we can formally define EOL with 


and RUL with 


EOL(t P ) = argmin T E OL(x(t),0{t)) = 1 , 

p 


RUL(t P ) = EOL(t P ) -t P . 

Due to process noise, measurement noise, and uncertainty in the future inputs of the system, we com- 
pute a probability distribution of the EOL or R.UL, i.e., at time tp, we compute p(EOL(tp)\yo-.t P ) and/or 
p(RUL(tp)\y 0 :t P ). 

II. B. Prognostics Architecture 

We adopt a model-based approach that utilizes detailed physics-based models of components and systems 
that include descriptions of how faults evolve in time. These models are parameterized by unknown and 
possibly time-varying wear parameters, 9{t). Therefore, our solution to the prognostics problem takes the 
perspective of joint state-parameter estimation. In discrete time k, we estimate x*, and 9 and use these 
estimates to predict EOL and R.UL when desired. 

We employ the prognostics architecture in Fig. 1. The system is provided with inputs u^- and provides 
measured outputs yu- Prognostics may begin at t = 0, with the damage estimation module determining 
estimates of the states and unknown parameters, represented as a probability distribution p(xfc, O^yo-.k)- 
The prediction module uses the joint state-parameter distribution, along with hypothesized future inputs, to 
compute EOL and RUL as probability distributions p(EOL^ p |y 0: fc p ) and p(RUL)- p |y 0: fc p ) at given prediction 
times kp. In parallel, a fault detection, isolation, and identification (FDII) module may be used to determine 
which damage mechanisms are active, represented as a fault set F. The damage estimation module may 
then use this result to limit the space of parameters that must be estimated. Alternatively, prognostics may 
begin only when diagnostics has completed. In this paper, we focus on the damage estimation and prediction 
modules, and begin prognostics at t = 0 without the aid of an FDII solution. 



Figure 1. Prognostics architecture. 


III. Valve Modeling 


Pneumatic valves are gas-actuated valves used in many domains. A normally-closed valve with a linear 
cylinder actuator is illustrated in Fig. 2. The valve is opened by filling the chamber below the piston with 
gas up to the supply pressure, and evacuating the chamber above the piston down to atmospheric pressure. 
The valve is closed by filling the chamber above the piston, and evacuating the chamber below the piston. 
The return spring ensures that the valve closes when pressure is lost. 

We develop a physics model of the valve based on mass and energy balances. The system state for the 
nominal model includes the position of the valve, x(t), the velocity of the valve, v(t), the mass of the gas in 
the volume above the piston, rrit(t), and the mass of the gas in the volume below the piston, m b (t): 


x(t) 


x(t) v(t) m t (t) 


1 T 


m b (t) 


The position when the valve is fully closed is defined as x = 0, and when fully open is defined as x = L s , 
where L s is the stroke length of the valve. 
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Figure 2. Pneumatic valve. 


The derivatives of the states are described by 


x(t) = 


v{t) a(t) f t (t) f b (t) 


1 T 


where a(t) is the valve acceleration, and ft(t') and f b (t) are the mass flows going into the top and bottom 
pneumatic ports, respectively. 

The inputs are considered to be 


u (t) 


1 T 


Pl(t ) Pr(t) U t (t) U b (t) , 


where Pi(t) and p r (t) are the fluid pressures on the left and right side of the plug, respectively, and u t (t) 
and u b (t) are the input pressures to the top and bottom pneumatic ports, which will alternate between the 
supply pressure and atmospheric pressure depending on the commanded valve position. 

The acceleration is defined by the combined mass of the piston and plug, m, and the sum of forces acting 
on the valve, which includes the forces from the pneumatic gas, (p b (t) — pt(t))A p , where p b (t) and pt(t) are 
the gas pressures on the bottom and the top of the piston, respectively, and A p is the surface area of the 
piston; the forces from the fluid flowing through the valve, (p r (t) — pi(t))A v , where A v is the area of the 
valve contacting the fluid; the weight of the moving parts of the valve, —mg, where g is the acceleration due 
to gravity; the spring force, —k(x(t) — x 0 ), where k is the spring constant and x 0 is the amount of spring 
compression when the valve is closed; friction, —rv(t), where r is the coefficient of kinetic friction, and the 
contact forces F c (t) at the boundaries of the valve motion, 


Fc(t) 


! k c (—x), if x < 0, 

0, if 0 < x < L s , 

—k c (x — L s ), iix>L s , 


where k c is the (large) spring constant associated with the flexible seals. 
The pressures p t (t) and p b (t ) are calculated as: 


Pt(t) 

Pb(t) 


m t (t)R g T 
Vt 0 + A p (L s — x(t)) 
m b (t)R g T 
Vb 0 + A p x(i) 


where we assume an isothermal process in which the (ideal) gas temperature is constant at T, R g is the gas 
constant for the pneumatic gas, and Vj 0 and V bo are the minimum gas volumes for the gas chambers above 
and below the piston, respectively. 

The gas flows are given by: 


ft(t) = fg(p t (t),Ut(t)) 

fb(t ) = f g (p b (t),u b (t)), 
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where f s defines gas flow through an orifice for choked and non-choked flow conditions.' Non-choked flow 
for pi > p 2 is given by f g , n c(pi,P 2 ) = 


CsAsPi^j Z pi gT ( 7 !i) ^(pi) (p') 

where 7 is the ratio of specific heats, Z is the gas compressibility factor, C s is the flow coefficient, and A s is 
the orifice area. Choked flow for pi > p 2, which occurs when P 1 /P 2 exceeds (^if) \ is given by 


fg,c{Pl,P2 ) — C s A s pi\j Z R gT ( 7 +lJ 
The overall gas flow equation is then given by 


fg(pi,P2) = { 


fg,nc(Pl 5 P 2 ) if Pi > P2 and ^ < (^±i) (T " 1) 

/g,c(Pi,P2) if Pi > P 2 and ^ 

-fg,nc(P2,Pi) if p 2 > Pi and gL < (T±i) 

[-fg,c(P2,Pl) if P 2 > Pi and > (2+1) OFnT 


iT 


We select our complete measurement vector as 

y (t) = x{t) Pt{t ) p b (t) f v (t) open(t) closed(t) 
where f v is the fluid flow through the valve: 

xCt') / 

fv{t) = —C v A v yJ^ I Pfi — p fr \sign(p f i -Pfr), 

C v is the (dimensionless) flow coefficient of the valve, p is the liquid density, and we assume a linear flow 
characteristic for the valve. The open(t ) and closed(t) signals are from discrete sensors which output 1 if 
the valve is in the fully opened or fully closed state: 

1, if x(t) > L s 

0, otherwise 

1, if x(t) < 0 
0, otherwise. 



Fig. 3 shows a nominal valve cycle. The valve is commanded to open at 0 s. The top pneumatic port 
opens to atmosphere and the bottom opens to the supply pressure. When the force on the underside of the 
piston is large enough to overcome the opposing forces, the valve begins to move upward as the pneumatic 
gas continues to flow into and out of the valve actuator. In a little over 1 s, the valve is completely open. 
The valve is commanded to close at 10 s. The bottom pneumatic port opens to atmosphere and the top 
opens to the supply pressure. When the force balance becomes negative, the valve starts to move downward, 
and completely closes within 1 s. The valve closes faster than it opens due to the return spring. 


III. A. Damage Modeling 

From valve documentation and historical maintenance records, we have identified friction damage, spring 
damage, internal valve leaks, and external valve leaks as faults that lead to EOL. The variables that describe 
the magnitude of damage augment the state vector x(f), and the parameters defining their progression form 
the unknown parameter vector 0(t). 

One damage mechanism present in valves is sliding wear, which results in an increase in friction. 8 We 
characterize friction damage as change in the friction coefficient, and model the damage progression in a 
form similar to the sliding wear equation: 2 

r{t) = w r \F f (t)v(t)\ 
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Valve Position 


Top Gas Pressure 


Open Indicator 


0.04 



5 10 15 

Time (s) 
Valve Flow 



5 10 15 

Time (s) 



53 0.5 



53 0.5 


5 10 15 

Time (s) 
Closed Indicator 


5 10 15 

Time (s) 


20 


20 


Figure 3. Nominal valve operation. 


where w r is the wear coefficient, and Ff(t) is the friction force. Friction damage causes the valve to both 
open and close more slowly. The valve is considered to be at EOL when it cannot open or close within 
required timing limits. We define r + as the maximum value of the friction parameter at which the valve still 
behaves within these limits. So, T E o i (x(t), 0(t)) = 1 if r(t) > r + . 

We assume a similar equation form for spring damage: 2 

k(t) = -w k \F a (t)v(t)\, 

where Wk is the spring wear coefficient and F s (t) is the spring force. The more the spring is used, the 
weaker it becomes, characterized by the change in the spring constant. Spring damage causes the valve to 
open faster and close slower. The spring must also remain strong enough to close the valve against system 
pressure when the actuating pressure is lost, and the loss of this function also defines EOL. We define k~ as 
the smallest value of k at which the valve will both close within the timing limit and fully close upon loss of 
supply pressure. So, Teol(x(<), 0(f)) = 1 also if k(t) < k~ . 

An internal leak in the valve can appear at the seal surrounding the piston as a result of sliding wear. 
The pneumatic gas is then able to flow between the volumes above and below the piston, decreasing the 
response time of the valve. We parameterize this leak by its equivalent orifice area, Aj(t), described by: 2 

Ai(t) = Wi\F f (t)v(t)\, 

where Wi is the wear coefficient. We define Af as the maximum internal leak area at which the valve opens 
and closes within the functional limits. So, TEOL{*-(t),0(t)) = 1 also if Ai(t) > Af . 

External leaks may also form at the actuator connections to the pneumatic gas supply, due to corrosion 
and other environmental factors. We assume the growth of the area of the leak holes, A e t (t ) and A e f,(f), is 
linear: 


Ae,t(f) ^e,t 

Ae,b(f) 

where the t and b subscripts denote a leak at the top and bottom pneumatic ports, respectively, and te e ,i and 
w e ^ are the wear coefficients. A leak at the top pneumatic port causes the valve to open more quickly and 
close more slowly, whereas a leak at the bottom pneumatic port causes the valve to open more slowly and 
open more quickly. We define the maximum leak hole areas at which the valve still opens and closes within 
the functional limits as A+ t and A+ b . So, TEOL(x(f), 0(f)) = 1 also if A eji (t) > A* >t or A e ^(t) > A* b . 
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Algorithm 1 SIR Filter 

Inputs: {(x l fc _ 1 ,e^_ 1 ),w^_ 1 }^| ; jUfc_i :fc ,y/ ! 
Outputs: {(xl,0 l k ),wl}^ =l 

for i = 1 to N do 

oi ~ P (Ok\oi-i) 

w\ <— p(yfc|xj., 0fe, u fc ) 

end for 

N 

i = 1 

for i = 1 to N do 
wt <- wl/W 

end for 

{(xfc,0fc),t4}A i <- Resample({(x| : ,6»l : ),Wfc}iI 1 ) 


IV. Damage Estimation 

In model-based prognostics, damage estimation is fundamentally a joint state-parameter estimation prob- 
lem, i.e., computation of p{~x. k , 0fc|yo : fc)- 3 Particle filters offer a general solution to this problem for non- 
linear systems with non-Gaussian noise terms, 9, 10 and have become a favored approach within prognos- 
tics. 2,3,6,11,12 As described in Section III, for the pneumatic valves, both the state equation f and the 
output equation h are highly nonlinear, and, further, h contains the discrete measurements open and closed, 
necessitating the use of particle filters. 

In particle filters, the state distribution is approximated by a set of discrete weighted samples, called 
particles. The particle approximation to the state distribution is given by 

with N denoting the number of particles. For particle i, x(. denotes the state estimates, 0 k denotes the 
parameter estimates, and w k denotes the weight. The posterior density is approximated by 

N 

p(x fc ,0 fe |y O :fc) « ^2w l k S { ^ t0 i k) (dx k de k ), 

i = 1 


with 0 ij(dx k dO k ) denoting the Dirac delta function located at (x(., 6).). 

We utilize the sampling importance resampling (SIR) particle filter, given as Algorithm 1. Each particle 
is propagated forward to time k by sampling new parameter and state values. The particle weight is assigned 
using y k . The weights are then normalized, and then the particles are resampled, e.g., using the systematic 
resampling algorithm. 13 

The parameters 0 k evolve by some unknown random process that is independent of the state x k . To 
perform parameter estimation within a particle filter framework, we assign a random walk evolution, i.e., for 
parameter 9, 9 k = 9 k - i+£,k-i, where £ k - 1 is a noise term. During the sampling step, particles are generated 
with parameter values that will be different from the current values of the parameters. The particles with 
parameter values closest to the true values should match the outputs better, and therefore be assigned 
higher weight. Resampling will cause more particles to be generated with similar values. The particle filter 
converges to the true values as the process is repeated over each step of the algorithm. 

The selected variance of the random walk noise determines both the rate of this convergence and the 
estimation performance once convergence is achieved. Therefore, this parameter should be tuned to obtain 
the best possible performance, but the optimal value is dependent on the value of the hidden wear parameter, 
which is unknown. We use the variance control method presented in Ref. 6, shown as Algorithm 2. In this 
approach, the variance of the hidden wear parameter estimate is controlled to a user-specified range by 
modifying the random walk noise variance. We assume that the £ values are tuned initially based on the 
maximum expected wear rates. The algorithm uses relative median absolute deviation (R.MAD) as the 
measure of spread, defined as: 


RMAD(X) = 100 


Median,; (| A, — Median., (X,)|) 
Median.,- (X,) 
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Algorithm 2 £ Adaptation 

Inputs: {(x^,0t),«4} iLi,Z k -! 

State: a 
Outputs: £ k 
if k = 0 then 
a <- 0 
end if 

for all j £ {1, 2, . . . , ng} do 

Vj - RMAD({^0')}£i) 

if a (j) = 0 and Vj < T then 

a (j) 1 

end if 

if a(j) = 0 then 

V* <- v* 0 

else 

Vj <— Wjoo 

end if 

CUD C, >-./*( l ! /’V') 

end for 


where A is a data set and X, is an element of that set. The adaptation scheme resembles a proportional 
control law, where the error between the actual R.MAD of a parameter 0(j), denoted as Vj in the algorithm, 
and the desired R.MAD value (e.g., 10%), denoted as v* in the algorithm, is normalized by Vj. The error is 
then multiplied by a factor P (e.g., 1 x 10 3 ), and the corresponding variance £(j) is increased or decreased 
by that percentage. There are two different setpoints. First, we allow for a convergence period, with setpoint 
v* 0 (e.g., 50%). Once Vj reaches T (e.g., 1.2i>* 0 ), we mark it using the a (j) flag, and begin to control it to a 
new setpoint (e.g., 10%). 


V. Prediction 


Prediction is initiated at a given time kp. Using the current state estimate, p(x^, 0k\yo:k) the goal is to 
compute p(EOL kp |yo : fc P ) and p{RULk P |yo : fc P ). The particle filter computes 

N 

p(xk P ,d kp \yo:kp ) tt^w l kp 8^ p ^i kp) (dx k pd0kp). 

i = 1 

We can approximate a prediction distribution n steps forward as 14 

N 

P(*k P +n? @kp-\-n |yO:fcp) « E u, *p 5 (*i P+ „.®ip+») (dXfc P+n^^fcp+n) • 

i= 1 

So, for a particle i propagated n steps forward without new data, we can simply take its weight as w l kp . 
Similarly, we can approximate the EOL as 


N 

p(EOL kp \y 0:kp ) « '£wi P 5 EOLip (dE0L kp ). 

i- 1 

To compute EOL, then, we propagate each particle forward to its own EOL and use that particle’s weight 
at kp for the weight of its EOL prediction. 

The pseudocode for the prediction procedure is given as Algorithm 3. 2 Each particle i is propagated 
forward until TporJ^l^ d‘ k ) evaluates to 1, at this point EOL has been reached for this particle. Prediction 
requires hypothesizing future inputs of the system, u^. For the valve, this problem is simplified, because 
each valve cycle corresponds to the same set of inputs, because the fluid pressures pl and pp can safely be 
assumed to be constant in our application domain, and, further, have an almost negligible effect on valve 
behavior because the forces they produce are very small compared to the other forces acting on the valve. 
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So, the future inputs for each valve cycle are deterministic. We can simply provide repeated valve cycles as 
input, and the prediction algorithm will determine EOL in the number of valve cycles. 

Fig. 4 shows results from the prediction of internal leak damage for N = 100. Initially, the particles have 
a very tight distribution of internal leak area values, A ; (f), but the distribution of the wear parameter, Wi, 
is relatively large. As a result, the individual trajectories for the different wear parameter values are easily 
distinguishable as EOL is approached. The different EOL values along with particle weights form an EOL 
distribution approximated by the probability mass function shown in the figure. The internal leak damage 
accrues step-wise as the valve moves, but over large time scales the damage progression appears linear, as 
shown in the figure. 

o.i - 

.£■ 

3 0.075 - 

ie EOL Probability Mass Function 

| 0.05 - 
£ 

0.025 - 


0-1 1 1 1 1 1 1 . 

50 60 70 80 90 100 110 120 

- 6 t (cycles) 



50 60 70 80 90 100 110 120 

t (cycles) 


Figure 4. Prediction of internal leak damage. 
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VI. Prognostics Interface Design 


Since prognostics is still a nascent field, there are no real examples of successful prognostics interfaces. 
Therefore, in this section, we try to establish some fundamental guidelines and principles for the design of 
prognostics interfaces, based on our experience developing a prototype interface for our case study. We first 
describe the information available from prognostics algorithms that may be presented to a user, then provide 
some design guidelines, followed by a description of the interface we designed. 

VI. A. Information Provided by Prognostics 

A prognostics interface must present the results of prognostics algorithms to a user and allow the user to 
interact with the algorithm to obtain additional results, view information in a variety of ways, and access 
long-term historical information. In contrast to most systems interfaces, e.g. diagnostics interfaces, in which 
information is viewed only over short time scales, such as the duration of a mission, prognostics interfaces 
cope with information over very large time scales, such as sequences of missions. Therefore it is important to 
distinguish these two time scales clearly. Further, it is more useful to an operator for results to be presented 
in an operational or missions context, e.g., how many missions a component will survive as opposed to how 
much time or how many cycles it will last. 

Prognostics algorithms are capable of presenting a wealth of information to the user. Foremost, algorithms 
generate predictions of EOL and RUL, i.e. , p(EOL(k p ) |yo : fc P ) and/or p(RUL(kp)\yo-.k P ) at given prediction 
times k p . Predictions may be computed automatically by the algorithm, but also on request from a user. 
As described in Section IV, model-based prognostics algorithms also compute supportive information, such 
as state and parameter estimates, p(x*,, 6^ |y 0: fe), which includes estimates of the amount of damage and the 
rates of damage progression. 

As shown in Algorithm 3, the state estimates are used directly to compute EOL. Using the same basic 
procedure, other predictions can be made as well, which may be useful to the user. Future state, parameter, 
and output trajectories can be computed and displayed to the user, which can show how the component 
is expected to perform in the future, and how damage is expected to progress. Note also that a number 
of health or damage indices may be computed from the state estimates, e.g., that a valve is at 50% health 
or its spring is at 75% health. Therefore, predictions of these health indices may also be obtained from 
the predictions of the state trajectories. This type of information is useful to the user because it helps in 
understanding the EOL and RUL predictions made by the algorithm. 

Once a component has failed, and ground truth is hence known, the outputs of the algorithm over the 
life of the component can be used to evaluate the performance of the algorithm. Various offline performance 
metrics may be used for this purpose, such as relative accuracy, convergence, prognostics horizon, and a-X 
performance. 15 This higher-level information is also beneficial to a system operator. 

VI. B. Design Guidelines for a Prognostics Interface 

With the wealth of information available from prognostics algorithms, a prognostics interface must avoid 
overwhelming the user with too much technical information. In the early life of a component, prognostics is 
not a concern, so the interface should indicate when damage has reached a threshold where EOL predictions 
become relevant so that the user knows when to start viewing the algorithm results. The most critical 
information, which should be displayed most prominently, includes overall component health and the life 
predictions. The predicted values should be shown over a large time window, so that users can see how the 
RUL and EOL predictions have evolved in time. Because it takes time for prognostics algorithms to converge, 
the interface should distinguish between predictions obtained before and after algorithm convergence, since 
the results before convergence are unreliable. 

State, parameter, and output estimates provide additional supportive information that helps the user 
to understand the current life predictions made by the algorithm. In particular, the evolution of damage 
states (e.g., spring rate or leakage area) over time demonstrate how the component is degrading, especially 
when shown against damage limits. It is equally important to show damage as it manifests in the functional 
performance of the component. For example, with the valves, EOL is largely defined by valve opening 
and closing times, so showing the progression of these values over time against the functional limits is also 
useful. Since damage often manifests as nonlinear changes in the system outputs over time (e.g., exponential 
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degradation), showing the predicted future output values helps the user understand the connection between 
functional performance and the life predictions with additional clarity. 

Since most prognostics algorithms provide results as probability distributions, there are multiple ways 
to view the results. Histograms or probability density functions can be shown, as well as box plots, or 
the life predictions at specified confidence levels. The choice of how to display the distribution depends 
on the type of reasoning that the user is performing. It is then also useful to show summary statistics 
for the various distributions, including measures of central tendency and statistical dispersion, as well as 
higher-order moments such as skew and kurtosis. 

Besides choosing how information should be displayed, the user may also interact directly with the 
algorithm by setting assumptions about the future inputs of the component. For example, a user may be 
interested only in worst-case performance, so will assume maximum loading on the component. Instead, a 
user may be interested in average performance over a distribution of possible future loading profiles. Each 
of these cases should be made available to the user. 

VI. C. Valve Prognostics Interface 

A prototype prognostics interface for the pneumatic valve case study was developed for demonstration 
purposes, and provides an implementation of the design concepts discussed in the previous sections. The 
interface is organized in a hierarchical structure that presents prognostics information at various levels, 
allowing the user to drill down into more detailed views when desired. 



Figure 5. Top-level systems interface with prognostic indicators. 


At the top level of the interface, shown in Fig. 5, is a system schematic integrated with sensor values. 
The only prognostics information communicated to the user at this level is the component health, shown 
using a simple color indicator with green indicating healthy, yellow indicating degraded health, and red 
indicating near-failure. This indicator changes continuously from green to yellow to red, in order to indicate 
the accumulation of component damage. When the indicator becomes yellow, this serves as a notification to 
the operator to investigate more closely. The next level of the interface, which is a component-centric view 
of prognostics results, allows this functionality. 

The lower-level prognostics interface, shown in Figs. 6 and 7, shows a summary of component health, the 
current life predictions, and supportive information. The number of fuelings the component has undergone, 
along with the overall health and RUL, are shown in text. The RUL is plotted on the left side, with the mean 
with confidence intervals shown in the figure. Here, multiple views are supported, including plots for the RUL 
at different confidence levels, the EOL mean with confidence intervals, the RUL probability distribution, the 
EOL probability distribution, and an a-X plot 15 that is available once the component has failed. Here, color 
is used to indicate the severity of the predictions, as the colors of the points are selected based on operational 
RUL thresholds. For example, the points are yellow when less than 15 remaining fuelings are predicted, as 
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Figure 6. Lower-level prognostics interface. 



Figure 7. Lower-level prognostics interface at component failure. 


at this point an operator might start to plan for future maintenance. The predictions up to the 10th fueling 
are shown in gray, as the algorithm has not yet converged until then in this example. 

The interface provides supportive information to the right of the RUL plot. At the top of this portion 
is shown the measured valve timing, which in this case is the only directly observable information for 
the selected component. The opening and closing time requirements are shown with red dashed lines. The 
expected future valve times are shown with the trajectories consisting of gray dots, with the mean, maximum, 
and minimum timing vales shown. These timing predictions help to explain why the algorithm predicts EOL 
at around 38 fuelings, because extrapolating the times linearly would yield a much larger EOL. Below the 
timing plots, the operator may select a specific damage mode, and view the estimated amount of damage 
and the estimated progression rate. Damage thresholds corresponding to EOL are shown with red dashed 
lines. 
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Figure 8. System report of prognostics information. 


A summary of prognostics information for multiple components, shown in Fig. 8, is also available. For 
each component, the overall health is shown, the RUL at the 95th percentile, and current opening and closing 
times, as well as the timing requirements. For each valve, the amount of damage for each damage mode is 
also shown. Damage modes with high damage percentages indicate where repairs will have to be performed. 
Such a report may be used to track the health of various components after each fueling operation. 

VII. Demonstration 

We apply the prognostics methodology to our case study, the transfer line chillldown valve of the Space 
Shuttle LH2 refueling system, using historical valve degradation data. Because the actual data consisted of 
only 7 valve cycles before maintenance was performed, we extend the data with simulated valve data, based 
on the identified bottom external leak damage mode with w e .b = 3.0 x 10“ ' m 2 /cycle, up to EOL at 38 
cycles, after which the valve fails to open within the 5 second limit. Fig. 9 shows the valve timing data used 
for the demonstration. Note that only the open and closed sensors are available here for prognosis, and we 
have established in previous work that prognostics may still be performed under this constraint. 2 




Figure 9. Real and simulated valve timing. 


We assumed that no FDII algorithm was present to limit the number of damage modes to estimate, and 
therefore damage estimation was started immediately, jointly estimating all the damage modes of the valve. 
Because of the large state space, we set N = 2000. The algorithm should converge to correct wear rates for 
each damage mode, and, in this case, identify the bottom external leak as the dominant damage mechanism. 
Estimation results are shown in Fig. 10. For the variance control algorithm, we used T = 60%, Vg = 50%, 
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Figure 10. Estimation of pneumatic valve damage modes. 


= 10%, and P = 2 x 10 -4 . The algorithm keeps the RMAD of the wear parameters controlled to within 
10%. The estimation algorithm correctly identifies the bottom external leak as the dominant damage mode, 
as shown by the convergence of its wear rate to the true value. The estimates for all other wear rates are 
driven to zero, and so the corresponding damage modes do not grow significantly. During the first 7 cycles 
(where measured data is available), the algorithm is still converging, and by the time this portion of the 
data has ended, the algorithm estimates the wear parameter for the bottom external leak within 30% of the 
true value, and the mean RUL prediction is within 25% of the true RUL. The predicted value of the RUL at 
that time point is less than the true value, resulting in a conservative estimate. By 8 cycles, the mean RUL 
prediction has improved to within 5% of the true RUL. At this point the algorithm has converged and the 
user may then begin to use the results in the decision-making process. 

Fig. 11 shows the RUL predictions at the 99% confidence level, i.e., the value at which 99% of the 
prediction distribution is greater than or equal to that value. The predictions at this confidence level always 
lie below the true RUL, therefore, they serve as conservative estimates upon which risk-averse decisions can 
be made. After algorithm convergence, the predictions are fairly smooth. The relative accuracy of the RUL, 
using the mean values, averages 91% over all predictions, and the RMAD of the RUL averages 9.88%. 
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Figure 11. RUL predictions at the 99% confidence level. Predictions before algorithm convergence are represented by 
unfilled circles, and predictions after convergence by filled circles. 


VIII. Conclusions 

In this paper, we developed a model-based prognostics approach, using particle filters for joint state- 
parameter estimation. The estimated damage state of the valve is propagated forward in time to predict 
EOL and RUL distributions. We adopted a pneumatic valve from the Space Shuttle refueling system as 
a case study. We developed a detailed physics-based model that includes damage progression, and applied 
the algorithmic approach to historical valve degradation data augmented with simulation data to EOL. A 
prototype user interface was also developed to communicate the prognostics results to a system operator. 
The results here demonstrated the effectiveness of a model-based approach for valve prognostics. 
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